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A am Mtlaata of jupit«r*« ploaotory flold Is obtsinsd fro* too 

VOyanor 1 obssnrstloas of too Jovisa •sdootosptam^o. An oxplioit mooI for tho 
MgaotoOiso ourrsat syatm is coablnsd wltli s sphoriosl Hsraonio aodol of tteo 
planstsry fiold with both sots of psrsMtors dotomiaod similtsaoousiy using s 
noa-linosr gonorsiisod invorso aotbodology. Iho rosulting aodol fits tho 
obsonrstioas ostroaoly wall throughout tho snslyais iatorval (r < 20 Jovian 
radii). Iho Jovian intomal fiold aoool obtsinod froa tho voysgor i oats is 
vory siailsr to tho ootopolo Pionoor 11 aodols. Iho boat fitting nsgnotodiso 
lios in tho oontrifugal oqustor, 2/3 of tho way botwoon tho rotational and 
aagnotie oquators, as appropriato tor oontrifugal loading of tho aagnotosphoro 
by a oold plaaaa. No statistically aignif leant ovidonco is found for secular 
change of tho oquivalont dipolo oatiaatod froa Pionoor 11 (197<i.9) and \fbyagor 
1 (1979.2) data. 


INTBODUCTIOM 

Iho Payagor i (VI) oncountor witn Jupitor in March 1979 was tho third of 
four such onoountors to provido dotollod in situ obsorvationa of tho Jovian 
aagnotosphoro. Iho low^latltudo approach of Pionoor 10 in Otcoabor 1973 to 
within 2.6 Jovian Radii planotocontrio distance (1 R. « 71323 «cb) provided tho 
first observations of tho Jovian aagnotodisc Chiith ot al.« 1974; Van Alien ot 
ai*» 1974) as well as tho first ostiaatos of Jupiter's internal aagnotie fioid 
based on in situ observations (Saith ot al., 1974). Tho high«latitudo , 
retrograde approach of tho Pionoor 11 (P11) spacecraft to within 1.6 R, in 
Oocoabor 1974 proved to bo tho aost favorable for tho estimation of Jupiter's 
internal fiold and led to spherical haraonio magnetic fiold aodols based on 
tho vector hoUum magnotomoter observations (Smith ot al., 1976) and ttie high 
fiold fluxgate magnotomotor obsorvtions (Acuna and Ness, 1976). 

I 

Preliminary attempts to obtain m internal field model from VI magneto- 

ootor observations (Ness ot al., I979) were frustrated by the large periapsis 

of VI relative to PIO and P11 (4.9 versus 2.8 and 1.6 R.) and the ubiquitious 
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pr«s«nG« of 0 largo aoalo equator iai currant systM aaaoclatod tilth the Jovian 
■agnotodlse. Thia diae-llko ayataa of oaatuard azlnuthal eurronts oxtends 
fraa Inaldo the orbit of lo (at 5.9 S.) outward to 50 I, and beyond 
(Gonnemey et al., 198I). Ihe 9.6* tilt of jupiter'a nagnetio dipole with 
respect to the rotation axis and the near-equatorial approach of Voyager 1 
resulted in the periodlo ianersion of the Itoyager spacecraft in the 
current-carrying region as it traversed the inner Jovian n^netop^ere. Ihus 
the traditional aethods of analysis of such data, utilizing <Mrthogonal 
spherical haraonic ^inctions to represent the aagnetic field, are not 
applicable since they require that the observations be obtained in a source 
free (current-free) region of ^>ace. Ibis is equivalent to the assuaption 
that the aagnetic field is derivable froa a scalar potential function. 

Connerney (1981) deaonstrated bow saall errors or unaodeled contributions 
to the observed aagnetic field (such as those due to local current systems) 
can lead to large errors in derived aagnetic field aodels. It is therefore 
essential to interpret the Voyager observations within the context of a model 
which is as representative of the physical situation as possible. Connerney 
(198I) introduced such a model for the analysis of aagnetic field observations 
at Ji^iter, incorporating explicitly the field contribution of large-scale 
external current systems in the Jovian magnetosphere, the observed field is 
modeled as the sum of two components. The planetary field is derivable from a 
scalar potential and represented by the usual spherical harmonic expansion. 

The external field, due to the distributed currents in Jupiter's magneto- 
sphere is derived from an appropriate vector potential. The parameters of 
both the Booel external current system and the model internal field are then 
determined simultaneously by inversion of the aagnetic field observations. 

From the Voyager 1 observations we are thus able to obtain an estimate of 
Jupiter's internal aagnetic field at epoch 1979.2 as well as a characteriza- 
tion of the aagnetodisc current system. We assume that the field of the 
external current system does not vary appreciably during the encounter period. 

Such a variation could possibly masquerade as a spatial variation of the 
field which would be reflected and not identified correctly in both sets of 
model parameters. 

METHODOUWY 
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Hm tMO •sMBtlAl (and aoval) f««tur«« of our anolysis of tho VMiygtr 
■•iMtio fiolO obsorvotlons aro (i) tOa usa of a aodal In which axtarnal 
eurraata ara aspiloitiy rapraaantaO and (2) tha uaa of ganaraliaad invarsa 
taehnlquaa to cOtain a aatlafaetory aolution. Both ara dasorlhaO In datall by 
Oonnarnay (1981) aa appliad to tha aralyais of Plonaar 11 ohaarvationa at 
Jupltar; wa wlU oitly ouUlna tha oathodology aa It appllaa to tha Voyagar 
ohaarvationa and aa it haa baan axtandad for tha praaant analyaia. 

lha aodal fiald ia rapraaantad aa tha aiat of an Intarnal fiaid B* 
darlvabia froa a aealar potantial and an axtarnal fiald b dua to aagnatodiac 
ourranta 


B a B* ♦ b* 

lha intarnal aagnatio fiald B* ia axpraaaad aa tha gradiant of a aealar 
potential function V. B* « whara 


V « a 


(a/r) i Pp“ (coa«) I g^^" coa(a^) ♦ h^^® ain(a^)) 
n«1 B«0 


r ia tha dlatanca to tha planat'a cantar, a ia tha planatary radiua, ^ and 
ara co-latitude and longituda* raapactivaly: the P^® are the aaaociated 
Legendre funetiona with Schaidt normalization, and the g^®, h^® (Sclialdt 
coefficienta) are the internal field parametera. 

Following Connerney et al. (1981) ue aaauaa that external aziauthal 
ourranta of tha aagnatodiac ara confined to an aziauthally aymaatric, 
planatocentrxc annular diac (Figure 1). The diac model parametera are the 
inner and outer edge radii and R^, the diac half-thicKneaa 0, and a acale 
conatant for the current denaity, which variea inveraaly with diatance from 
Jupiter. Two additional free parametera and apaoify tne orientation of 
tha current diac with reapect to Jupiter (Syatea 111 1965): the normal to tne 
current diac makea an angle with Jupiter 'a rotation axia and liea in tha ^ 
(weat) longitudinal meridian. The external field b due to theae currenta ia 
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ooaput«4 RiMwrioally 


Uith th« addition of ttaa aagnacodiao ourront ayataa* tha aooal nagnotie 
fiald is no longar linaar in tha nodal paranatars and wa aust uaa an itarativa 
invarsion taehniqua daacribad by Connarnay (19d1)« appropriataly aodifiad to 
aoooaaodata tha paranatars of tha eurrant disc as fraa paraaatars. Iha 
linaarizad systaa to ba solvad ac aaoh itaration is 

y « Ax 

uhara y is a ooliaui vactor of tha aodal rasiduals (ttta obsarvad ainus aodalad 
fiald), X is a coltaui vaotor oonsisting of tha paraaatar corractions raquirad 
to bring tha aodal into olosar agraaaant with tha data, and tha aatrix A is a 
aatrix of partial darivativas of tha aodal fiald with raapact to tha model 
paraaatars. Tha vactor y is of length M, x is of length M, and A is an N by M 
aatrix, where N is the nuaber of (coaponent) magnetic field observations, and 
H is tha number of fl*ee paraaeters 

M « (n^ ♦ 1)^ - 1 ♦ 6 

associated with tha internal field expansion to order n and the 6 

max 

paraaatars of the current disc. Inclusion of the disc parameters R R,. 0, 

and requires a transformation (sealing) of the parameter vector (not 
discussed by (k>nnarney (1981); see, e.g., Lawson and Hanson (1974)): 


where is tha axpaotad standard deviation of the parameter correction x^. 

We adopt a relative scaling of * 1 for the internal field parameter 

corrections (4g ", Ah ") and suitably chosen for the corrections to the 
n n 1 

disc paraaeters such that in program units the components of the solution 
vector X are approximately equal in magnitude. Additionally tne partial 
derivatives of the model with respect to the parameters of the mooel current 
disc (rightmost 6 columns of the A aatrix) must be computed numerically, as is 
the modal field. 
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With ttt«M modifieationst tb« mtbod of oonstruotlog solutions outlinsd 
in Connornty (1981) is usod. Wo ehooso on intornal sphoriool harmonic 
expansion of order n^^ s 3 to facilitate direct and meaningful comparisons of 
VI internal field models with previous models (of internal order 3} obtained 
from PII observations. Our model thus has 21 free parameters, not all of 
which will be determined from the available VI observations; Uie 
interpretation of such insufficient data requires the construction of partial 
solutions. The singular value decomposition of Lanosos (Lanezos. 196I) is 
used to reformulate the problem in terms of independent parameter vectors 
(eigenparameters) which are linear combinations of the original model 
parameters. A solution is constructed by summation over a subset of the 
eigenparameters. starting with a few well determined parameter vectors and 
successively increasing the number of eigenparameters in the solution. As 
additional eigenparameters are added, more of the original model parameters 
are resolved; eventually reaching a point where the remaining paramete** 
vectors are so poorly determined (due to the limited observations) that 
inclusion seriously degrades the solution. The partial solution eonstructeo 
xn this way represents the best available estimate of the solution, ano it is 
understood that the remaining parameter vectors (not used in the solution) are 
undetermined . 

RESULTS 

Our Voyager 1 data set selected for inversion consists of 500 vector 
observations of the magnetic field taken every 6 minutes curing the interval 
from hour 16 day 63 to hour 8 day 65* during whicn the radial distance of 
Voyager 1 from Jupiter ranged from 20 R, to the close approach (hour 12. day 
64} distance of 4.9 R,. Each observation is a 48 second average of vector 
observations obtained every 60 msec with an estimated accuracy of 0.2 nT ^0.1} 
of full scale (Ness et al.. 1979). The more distant observations tnat are 
less sensitive to the internal field parameters are included because they 
greatly improve the resolution of the current disc parameters. Improved 
resolution of the current disc parameters leads to a more confident separation 
of internal and external (local) fields and ultimately an improved internal 
field model. More traditional analyses of such data that do not explicitly 
include models of the external current system (e.g., Acuna ano Ness (1976); 
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AolUM «t (1981); aaith «t «1., (1976)) liait th« obsarvations us«d to r < 
8 or 10 Rj In an attaapt to ainiai :• t)M iapaot of auoh axtornal fiold 
oontrlbutlona . 

Iho Itarativo invaralon tachnlqua raquiroa an initial paraaatar a«t about 
Mbieb tba problaa ia eonaidarad to ba auffioiantly linaar locally that 
auoeaaaiva applioationa of tba linaar genaralizad invaraa tacbniquaa raault in 
a oonvargaot aolution. Ua aalaetad aa initial aodala aiapla til tad dipola 
internal fialda (n s 1 taraa) ooaoinad with tba VI nodal magnatoapbara currant 
diao paranatara (Gonnarnay at al. (1981))* Savaral dipolar intarnal field 
■odala, intentionally diaplacad from the lowaat oroar dipola obtained from P11 
obaarvationa • ware uaad to damonatrata tbat tba final aolution did not depend 
on tba initial nodal. A typical initial intarnal field nodal ia cbaractarized 
by tba paranatar aat s M.O G, g^^ s .3 G. b^^ « .3 G« g^^* * 0 and b^* s 0 
for all n > 1, corraaponding to a ainpla dipole tilted by 6* towarda a Syatem 
III longitude of 225*. By conpariaon, the Pli dipoles are typically tilted by 
'*‘10* towards of ^ 200'*. 

After 12 iterations we obtain the solution listed in Table 1. No entry 
is made in the table for the nodal paraneters that are unresolved (defined as 
parameters with corresponding resolution matrix elements of < 0.95; (see, 
a.g., Connerney (1981); Wiggins (1972); or Jackson (I972)). For comparison we 
list in Table 1 several models based on the Pioneer 11 observations; in 
general an excellent agreenent is found between the VI parameters and chose of 
the P11 based models. The unweighted RMS of the VI model residuals tnroughout 
the entire data interval of R < 20 R, is 7.8 nT. The quantity < listed with 
the VI model internal field parameters corresponds to an estimated 2° error 
assuming uncorrelated errors; the true estimated errors are certainly greater 
than c since the errors are in fact correlated. But the quantity * is 
expected to give some indication of the relative errors among tne parameters 
and are presented for that purpose. 

The resulting parameters of tne current disc are not unlike those quoted 
by Connerney et al. although that model was not an optimal fit to the 
observations. The most interesting difference between the two is that tue VI 
optimal fit yields a current disc which is not coincident with the magnetic 
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•qiMtor (tilted by 9.6* from the rotetlon «x«b) but inatbad ia tutad by 6.5*« 
approxiaataly 2/3 of tba uay batvwan tha rotation axia and tha aagnatio 
equator. 

in illuatration of how wall tha nodal fita tha Voyager l obaarvationa la 
ahown in Figira 2. lha perturbation field tg ia tha difference between the 
obaerved nagnetie field at any poaition and the field of internal origin aa 
obtained fron tha nodal fit. lha daahad line ia tha field of tha modal 

Qurrant diac; tha difference between the obaarvationa and the daahad line 

rapraaanta tna raaiduala, i.e«, tha nodal niafit. Thia rapreaantation ia 
ohoaan to anphaaixa tha relative nagnituda of tha field due to tha local 
nagnatodiac curranta, and it facilitatea an interpretation of Uia ranaining 

nodal raaiduala. lha reaiduala are vary anall for r < 10 noat of tha 7.8 

nT BHS raaiduala appeara at larger radial diatancea where tha field of tha 
external curranta ia a large fraction of the total field, lha raaiduala for r 
< 10 Rj are an exceedingly anall fraction of the total field, which growa co 
3330 nT at cloaa approach (4.9 Rj). lha vary large and localized feature 
evident at day 64 hour 15 ia the aignature of the intenae currant ayatem 
generated by tha interaction of the Jovian nagnetoaphere with tha aatellite lo 
(Kaaa at al., 1979; Acuna et al.. 198I). 

DISCUSSION 

lha modal of Jupiter 'a internal magnetic field at epoch 1979.2 ootained 

from the Voyager 1 obaarvationa beara a very cloae relationahip to the epocn 

1974.9 mooela obtained from Pioneer 11 obaarvationa. In general, the 

paranetera (e.g., g^\ h^^) tnat are expected to be relatively well determined 

are indeed the noat oonaiatent. lha cloae correapondence between tne dipole 

terma (n « 1) of the Voyager 1 model and the (Pioneer 11} 0^ model ia 

particularly atriking. Ihe VI and 0^ dipole tema g^^. g^^ and h^^ differ by 

only 0.24, 0.61 and 1.2S, reapectively . ^eae differencea are much amaller 

than the eatimated parameter uncertaintiea . It ia Intereating to note that 

the voyager 1 paracetera which do not agree aa well with the Pioneer li moUeia 
0 3 

(g^ and g^**) are alao the aubject of aome diaagreement among the Pioneer 
nodela obtained fron the two magnetic field experimenta onboard Pioneer 11. 

In compariaon, the preliminary eatimatea of Jupiter 'a magnetic field obtained 
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from m eonv«ntioo«l •ptatrleal iwraonlc analyais of tho Voyagor i obaorvatlona 
(Naaa ot aX., 1979) ylaldad oatiaatoa of jupitor*a dlpolo tilt ranging from 
9»6* to 13.?* towarOa longitudoa of 189* to I94*; oatlaatoa of tho dlpolo 
■agnltudo ranged froo 3,76 G»Rj^ to 4.09 G-Rj^. PIO aodola (Salth ot al., 
1976) boar ouob loaa roaonblanoo to olthor tho P11 or VI oodolt particularly 
In tho higher order ooofflelonta. It la clear that tho Pioneer 10 obaorva- 
tlons wore alao heavily Influenced by the magnetodiao currenta (Connerney et 
al.» 1981). 

Aaong the unreaolved VI Bodel paraaetera» the inner current aheet 

edge, and g^^ will alaoat certainly never be obtainable froa the Voyager 1 

obaorvatlona alone. Iheae two paranetera are heavily repreaentec In the aoat 

poorly deteralned eigenvector (21)* which la <^2 ordera of aagnltude aorc 

poorly deterained than any included in the Voyager l aolution. That la* 

within the context of the choaen pnyaloal aodel, the obaorvatlona are alaply 

Inaufficlent to deteralne the valuea of theae paraaetera* aa a conaequence of 

the apatlal dlatrlbutlon of the obaorvatlona. The difficulty of deteralnlng 

the Inner edge of the current aheet (R^) froa Voyager 1 ooaervaclona waa 

deduced intulclvely by Connerney et al. (1981) and la conflraed by the 

generalized inverae analyala. Independent obaorvatlona* however, auggeat that 

R la indeed cloae to 5 R,: the near axla external field deduced froa P11 

obaorvatlona (Salth et al., 1976) la very cloae to that expected of the 

current dlac of R 5 R,* The reaalnlng unreaolved paraaetera g ^ and h ^ 
o j 3 j 

are aaaoclated with eigenvectora (18 and 19) that are not aa dlacouraglng aa 
the aoat poorly determined eigenvector. It la conceivable that further 
analyala may provide at leaat aome Information about these parameters. 

The most Interesting result of tne optimal fit to the Voyager 1 
observations with respect to the current disc in Jupiter’s Inner magnetospnere 
la its orientation. The Voyager 1 observations at radial distances of less 
than 20 R^ are beat fit by a current disc not in the aagnetlc equator, as 
argued by Connerney et al. (1981)* Goertz et al. (1976) and Goertz (1976, 
1979)* but rather a current disc residing in a plane txlted only <’'2/3 of the 
way towards tne aagnetlc equator. Connerney et al. (198I) noted that the 
distant VI observations were insensitive to the disc orientation parameter ^ 
and argued aa did Goertz (1976, 1979) » that the current sheet resides close to 
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th« MgMtio tquator on th« basis of PIO obsarvationa 


nrior to any of tba Jupltar anoountars, dadhlll (1967) pradiotaO that 
tha oantrifugal foroa dua to Jupltar* a rapid rotation would confina a plasma 
to a diao-ahapad region in a plana tilted by 7* to Jupiter's equator. Hill at 
al. (1974) referred to that plana as tha 'centrifugal symmetry surface* to 
vfhlch cold plasma would be confined (sea also (joartz, 1976). For a hot 
plasma, tha pressure gradient and magnetic mirror forces dominate tha 
centrifugal forces (e.g., (toertz, 1976) and the plasma would reside in the 
magnetic equator. Thus it would appear that, within the context of our model, 
these results require the current in Jupiter's (inner) magnetosphere to be 
carried by 'cold* and not 'hot* plasma. However, it is precisely at this 
level of interpretation that the limitations of our current disc model arise. 
In particular, the disc thickness is assumed constant in radial distance, and 
the model azimuthal current is distributed uniformly in z. While such a model 
is capable of fitting the observations exceedingly well, it is possible that 
an equally good fit can be obtaine. with an alternate model. It may be 
possible, for example, to adjust the distribution of current in z within the 
disc and the disc orientation to obtain a model field similar to that 
illustrated in Figure 2 but with a disc oriented in the magnetic equator. 

Until the physical validity of our present current disc model can be 
ascertained by a self ..consistent treatment of the Jovian plasma and magnetic 
field, we regard the inferred orientation of the current disc as tentative. 

CONCLUSIONS 

The kind of model applied herein to the Voyager 1 observations, in wnich 
an internal spherical harmonic expansion is combined with an explicit model of 
the field due to external current systems is regarded as essential to 
understanding aiid integrating the magnetic field observations of each of the 
Jovian encounters, ^deed, the success of the model used is a very 
encouraging indication of tne extent of present knowledge of Jupiter's 
magnetic field (and external current system). We obtain from the Voyager 1 
data a Jovian internal field model for epoch 1979*2 that is independent of the 
previous Pioneer 11 observations and quite consistent with the epoch 1974.9 
Pioneer 11 models. 


10 



Ih« Voy«g«r 1 internal field model deduced Here atiould provide a basia 
for a rational diaouaaion of a poaaible aeoular variation of Jupiter 'a 
internal field. The atr iking reaemblance between the VI and 0,^ modela 
auggeats that Jupiter's internal magnetio field has not changed between the 
Pioneer li enoounter in Deoember 1974 and the Voyager 1 encounter in March 
1979* for example, equivalent tilted, centered dipoles ot the VI and 0^ 
models differ by only 0.2S1 in magnitude, 0.04* in tilt and 0 . 1 * in longitude 
and these are much smaller than the estimated parameter uncertainties. Thus 
we find no statistically significant evidence for any secular change. Our 
goal was to provide the beat independent estimate of Jupiter's internal field; 
a combined fit to various of the data seta available may yield an impiuved 
internal field model provided the observations themselves (Voyagers 1 and 2; 
Pioneers 10 and 11 } can be sensibly integrated. 
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.483 

(.018) 

.513 

(.048) 

.502 

.537 


7 

“a’ 

-.294 

(.058) 

-.469 

(.037) 

-.498 

-.531 

NA 

8 

-a^ 

.107 

(.018) 

.088 

( .037) 

.119 

.074 


9 

< 

— 

( — ) 

-.233 

(.060) 

-.111 

-.113 


10 

»3’ 

— 

( — ) 

-.076 

(.083) 

-.316 

-.585 


11 


.263 

(.114) 

.168 

(.080) 

.220 

.283 


12 

«3^ 

-.069 

(.054) 

-.231 

(.082) 

-.250 

.067 


13 

•’3’ 

— 

( — ) 

-.580 

(.104) 

-.476 

-.423 


14 

-3* 

.695 

( . 108) 

.487 

(.108) 

.380 

.120 


15 

-3^ 

-.247 

(.054) 

-.294 

(.068) 

-.228 

-.171 


16 

■0 

— 






5 1 
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17 

«1 

56 




50 Rj 

18 

0 

3.1 

N/A 

M/A 

N/A 

2.5 Rj 

19 

I 

185 




225 S/F 

20 

• 

0 

6.5 




9.6 DEG 

21 

♦0 

206. 




202. DEG 

1965 

Systen 

Ill ♦ 
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Schaictt Doraalized spherical haraonio coefficients. Gauss 

A. icuM and Ness (I 976 ): Rotated to 1965 Systen III 

B. Connerney (1981) 

C. Siiith et al. (1976: Rotated to 1965 Systeo III 

0. Connerney et al. (I 98 I) 
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FXGUKE 2 


A orosA-Motlon of tho Mdtl eurrtnt Nloe. 

OMporlson of Mdoitd porturbotlon Mcootlo flolo (4«sh«<l) with 
that ^oorvod for Voyogor i (aphorleol eoordiontoa art uaoo) * In 
thla proaontatlon tho aodol iatornal flold tiaa boon aubtraotoo fron 
tbo obaorvatlona : tho total flold at eloaoat approach (C.A.) la 
3330 nl* 
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FIGURE 1 




VOYAGER 1 JUPITER PERTURBATION FIELD 
^■OBSERVED- INTERNAL FIELD MODEL 






